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Abstract. A sample of metal-poor subgiants has been observed with the UVES spectrograph at the Very Large 
Telescope and abundances of Li and Be have been determined. Typical signal-to-noise per spectral bin values for 
the co-added spectra are of the order of 500 for the Lii line (670.78 nm) and 100 for the Bell doublet lines (313.04 
nm). The spectral analysis of the observations was carried out using the Uppsala suite of codes and MARCS 
(ID-LTE) model atmospheres with stellar parameters from photometry, parallaxes, isochrones and Fe ll lines. 
Abundance estimates of the light elements were corrected for departures from local thermodynamic equilibrium 
in the line formation. Effective temperatures and Li abundances seem to be correlated and Be abundances correlate 
with [O/H]. Standard models predict Li and Be abundances approximately one order of magnitude lower than 
main-sequence values which is in general agreement with the observations. On average, our observed depletions 
seem to be 0.1 dex smaller and between 0.2 and 0.4 dex larger (depending on which reference is taken) than 
those predicted for Li and Be, respectively. This is not surprising since the initial Li abundance, as derived from 
main-sequence stars on the Spite plateau, may be systematically in error by 0.1 dex or more, and uncertainties 
in the spectrum normalisation and continuum drawing may affect our Be abundances systematically. 
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1. Introduction 

Stellar abundances of the very light elements lithium, 
beryllium, and boron (Li, Be, and B) have been studied 
extensively during the recent decades, mainly because of 
the great role they play in constraining Big Bang nucleo¬ 
synthesis theories as well as models of galactic chemical 
evolution and stellar evolution. 

Stellar evolution predicts depletion of the light ele¬ 
ments in stellar atmospheres under certain physical con¬ 
ditions. Stars with their surface convection base reaching 
depths where there has been some depletion should re¬ 
flect such depletion in their stellar surfaces. The convec¬ 
tive motions can bring up material from the interior to the 
surface, i.e. gas poor in light elements is mixed with gas 
with the initial composition, leading to dilution, this is ex¬ 
pected to happen in subgiant stars which are the type of 
stars under study in this work. If the temperature at the 
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base of the convection zone is high enough, light elements 
brought down from the surface can be destroyed. 

With this background in mind, we have looked at 
lithium and beryllium abundances in a sample of metal- 
poor subgiant stars, from a new set of UVES observations. 
The analysis of the Be features, in crowded spectral re¬ 
gions where blending may be a serious problem, requires 
detailed spectrum synthesis of high quality. 

Standard models, e.g. lPelivannis et al.1 lll99fl(l . include 
only convection as the mixing mechanism. According to 
these models, depletion occurs at low stellar mass and 
low effective temperatures. These stars have convection 
zones deep enough to reach inner layers where the light 
elements have been depleted. The models predict burning 
of the light isotopes during the pre-main sequence (PMS) 
phase and dilution in the subgiant evolutionary phase. The 
fact that the Sun has a photospheric logarith mic Li abun¬ 
dance of A(Li) = 1.05 llAsplund et al.ll200^ which is sig¬ 
nificantly lower than the values measured in meteoritic 
matter (which is thought to be representative of the gas 
at the time our solar system formed) indicates that also 
main-sequence (MS) stars undergo some mixing process 
(other than simple convection). 
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O n the main-seque nce, diffusion I Michandl 
|l986LlRichar^ |200^, rotation i nduced mixing 

Til 99ft. mass l oss llSchramm et al.l 
Il990t I ^udai^^^Sarbonn^ [l995|) and gravity waves 

llOarcia Lopez fc 'sDmitni99iriTalon fc Charbonne]ll2004^ 

have been proposed to explain the light element patterns 
observed in MS stars. Light elements burn at low, but 
different, burning temperatures. The most fragile of them 
is Li (with a burning temperature Tb ~ 2.5 x 10® K), fol¬ 
lowed by Be (Tb 3.0 x 10® K) and B (Tb ~ 5.0 x 10® K). 
Therefore, if at least the abundances of two of the light 
isotopes are known in a star, their abundances may be 
used to further constrain which type of internal mixing 
is active in the stellar external layers. Models with only 
mass loss fail to predict a simultaneous Li and Be deple¬ 
tion, i.e. lithium is expected to be completely depleted 
before beryllium shows any sign of depletion. However 
observations show just the opposite, i.e. some degree of 
Be depletion in stars in which Li can still be detected. 
One such example (there are several other cases) are the 
Hyades F stars, which show depletion in both Li a nd Be 
jBoesgaard fc Tripicc^IlQSdt IPdivannis et ai]|l998ll. 

As for the case of main-sequence stars, the knowledge 
of two or more light isotopes in subgiants can help us 
to better understand the mixing processes occurring in 
the atmospheres of these stars. In this work, we have ex¬ 
ploited that possibility by determining both Li and Be 
abundances and comparing the observed abundance de¬ 
pletions with predictions based on standard models. 

This comparison requires the assumption of a value for 
the initial abundances so that the magnitude of the ob¬ 
served depletion can be quantified. About 70% of the ’’Li 
observ ed in the early Galaxy is synthesised during the Big 
Bang llWalker et alJ ll99ll Il993t lOlive et al.l l2o'o^. Some 


extra contributions by a — a reactions during late stellar 
evolutionary stages, like (post-)AGB stars and novae, con¬ 
tribute at later epochs to the enrichment of the Galaxy. 
These sources influence the galactic evolution of lithium at 
higher metallicities. Therefore, ’Li abundances in the old¬ 
est, most metal-poor stars of the Galactic halo are of great 
importance from a cosmological point of view, providing 
useful constraints on the co smological baryon-to- photon 
ratio T], In the early 1980 ’s, ISoite fc Soitei lll982h found 
a constant logarithmic lithium abundance, ^(Li) = 2.1 
in a small sample of warm metal-poor main-sequence 
stars spanning a wide range of (low) metallicities, the so 
called Spite lithium plateau. Although the plateau has 
now been confirmed by large data samples, the interpre¬ 
tation of its constant value (whether Li is primordial or 
not) still remains an open question, as well as if it is tilted 
with respect to metallicitv and effective temperature (cf 
iGharbonnel fc Primal l200^ for a recent review). If this 
plateau value represents the primordial abundance of Li, 
then one needs to explain why in the Galaxy today its 
abundance is 3.25 llAsplund et alJlio?)^ . 


Abundances are denoted by A(X) = log (Ax/A h) + 12 


For the other light elements, iReeves et al 1 lll97(lll and 
iMeneguzzi et al ] lll97lll proposed cosmic-ray spallation re¬ 
actions as the main (if not unique, as in the case of 
beryllium) source of ®Li, ®Be and production. In 

this scenario, light elements are produced via collisions 
of G and O target nuclei with protons and a particles in 
the interstellar medium. This predicts a quadratic depen¬ 
dence of the light isotope abundance on th e target nuclei, 
e.g. ox ygen. Howeve r, in the early 1990’s , iGilmore et alJ 
1 I 1992 I) ‘ (for Be) and IPuncan et al ] lll992h (for B) found 
a linear slope, which has triggered the exploration of 
several variations of the standard cosmic-ray spallation 
theory. The observ ed linear trend (further confirmed by 
larger data samples, [Moia^^^O 

lll999ll 'l seem to imply a primary origin of Be and B in¬ 
stead of a secondary one. One possibility to achieve this is 
if G and O nuclei were freshly synthesised by Type H su¬ 
pernovae (SNell) and accelerated by the outwards shocks 
together with the a particles and protons, and different 
scenario s more or less built on this idea have been ex¬ 
ploited ([Duncan et alj|l992| [Feltzin g fc Gustafssonlll994t 
iGasse et alJll99,5t IParizot fc Drurvlll999ll . Knowledge of 

the correlation between Be and O stellar abundances in 
larger samples of stars will help us to discern between 
these different scenarios (it should be noted that many 
more models have been developed in the last few years, 
and that those mentioned above represent some main 
ideas, which have subsequently been developed with dif¬ 
ferent degrees of sophistication). 

On the observational side, in general, the abundances 
of light elements in stars are so small that the only lines 
used for abundance determination purposes are resonance 
lines. The observations of these lines as well as their mod¬ 
elling have been quite a challenge. Most studies have so far 
concentrated on lithium; its resonant doublet at 670.78 nm 
makes it easily accessible with ground based telescopes. 
Beryllium, instead, having its main resonant doublet at 
313.0 nm, requires high resolution and highly efficient 
spectrographs in the near-UV. Boron, the third light ele¬ 
ment, is even more challenging as it may be studied only 
from space (such as with the Hubble Space Telescope), as 
its main resonant doublet falls at 250.0 nm. Thanks to 
the very efficient high-resolution spectrographs currently 
available, we note that observations of ®Li in Galactic stars 
have become more accurate and affordable in term s of 
telescope time llNissen et al ] l2000t lAoki et al ] |2004 and 
references therein). 

In the following sections we describe the observations 
as well as the stellar parameters and stellar abundances 
determinations, which are then discussed in the last sec¬ 
tion in terms of dilution in subgiants. 

2. Observations and data reduction 

Our sample of metal-poor subgiants was observed in 
service mode between October 2000 and March 2001, 
with the Ultraviole t Visual Echelle Spectrograph (UVES, 
iDekker et all [20^, mounted on the Nasmyth B plat- 
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form of the 8 m Kueyen unit of the ESO Very Large 
Telescope (VLT). All observations were taken in dichroic 
mode, i.e. using both the red and the blue arms of the 
spectrograph simultaneously, thus covering a wide range 
in wavelength. The instrument set-ups chosen were the 
Dichroic#! (Blue346 combined with Red580, where 346 
and 580 denote the central wavelengths in nm of the two 
settings), and Dichroic #2 (Blue346 plus Red860). In this 
way we obtained a spectral coverage from 300 to 380 nm 
in the blue, and from 480 to 1000 nm in the red (with the 
only exception of the narrow gaps between the two CCDs 
used in the UVES Red Arm). 

The observations had been proposed and planned for a 
detailed study of oxygen abundances using three different 
abundan ce indicators ([Oi], Oi, an d the OH lines in the 
near-UV, ICarcia Perez et al ]|2^ . but the large spectral 
coverage also allowed us to study other elements like the 
light elements, lithium and beryllium. The need for de¬ 
tecting the three different oxygen indicators in the same 
object restricted this study to a relatively narrow range of 
surface gravities and effective temperatures (see Tabled), 
which corresponds to stars that have already left the MS 
and started evolving (ascending) along the subgiant (gi¬ 
ant) branch. Moreover, accurate oxygen abundances de¬ 
terminations require accurate surface gravity determina¬ 
tions. Hence, our stars were chosen based on the availabil¬ 
ity of Hipparcos parallaxes from which surface gravities 
were estimated. The selected stars ended up being not 
very distant objects (the parallaxes of most of them in¬ 
dicate distances smaller than 400kpc) which reduces the 
effective temperature uncertainties due to reddening. 

Lithium is easily measurable from its resonant doublet 
at 670.78 nm, whereas beryllium is more challenging as 
its main resonant doublet falls at 313.0 nm, i.e. in a very 
crowded region and very close to the atmospheric cut-off. 
In order to maximise the signal-to-noise ratio of the near- 
UV spectra, a 2 x 2 binning was selected when reading out 
the EEV chip of the UVES blue arm. This choice did not 
affect the resolution of our final spectra since we chose 
a slit width of the nominal resolving power of which 
[R = 45 000) implies some oversampling. In the red arm, 
we selected a slit width of 0.7"corresponding to a nominal 
resolving power of i? = 55 000, thus resolving the blending 
of the Li line by an Fei line at 670.74 nm (in spectra 
for which such blending could be a problem, i.e. for the 
more metal-rich stars). Typical signal-to-noise {SjN') per 
spectral bin values are of the order of 100 and 500 at the 
Be and Li lines, respectively, in the final co-added spectra. 

Data reduction followed standard reduction proce¬ 
dures, i.e. bias subtraction, division by a normalised flat- 
held, extraction of the echelle orders and their wavelength 
calibration. After a careful inspection of the blue spectra, 
we decided that the quality achieved by the UVES data 
reduction pipeline processing was satisfactory for our pur¬ 
poses. However, for the red spectra, we decided to man¬ 
ually re-reduce all the spectra (using IRAF) in order to 
improve (i.e. minimise) the presence of fringes and ripples 
in these very high SjN spectra. We then normalised all 




Fig. 1. Spectra of the Be (top) and Li (bottom) regions 
in one star of our sample, HD 218857. The positions of the 
Be II doublet and of the Li i line are indicated by vertical 
bars in both panels. 


the spectra with a spline function. This step is straightfor¬ 
ward in the lithium region (where the stellar continuum is 
easily identifiable), but the severe spectral-line crowding 
in the near-UV region does not allow such an easy iden¬ 
tification. Instead, the near-UV spectra were normalised 
by using synthetic spectra that were computed for each 
individual star. 

In spite of the weakness of the observed lines, the very 
high quality of our UVES spectra allowed us to detect both 
the Li and the Be lines in almost all stars. Fig.^illustrates 
the spectral quality of one of our targets, HD 218857 (U = 
8.967). 

3. Stellar parameters and model atmospheres 

3.1. Stellar parameters 

Since the determination of the stellar parameters for 
our sample of stars is presented in a parallel paper 
I Carcfa, Perez et uere we will only summarise 

how they were derived. Table ^ lists the values used for 
the effective temperature (Tefi), the logarithmic surface 
gravity (logg), the metallicity ([Fe/H]) and the microtur¬ 
bulence. 

The process of determining stellar parameters is clearly 
an iterative process. In short, T^s was estim ated from the 
b — y )- and (U — K)-Teff calibrations of I Alonso et alJ 
1996li . based on the Infrared Flux Method (IRFM). 
Stromgren (uvbyP) and V, K photometry were taken 
from Schuster & Nissen (1988, and private communica¬ 
tion) ) and the Two Micron All Sky Survey (2MASS) 
database, respectively. Colours and magnit udes were de- 
redden ed with values estimated from the iHakkila et alJ 
(1993) galactic maps of interstellar reddening. 
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The surface gravity was derived from Hipparcos par - 
allaxes following the prescription in iNissen et alJ lll997ll . 
Parallaxes are available for all programme stars and, fur¬ 
thermore, their uncertainties are all very small, except 
for two stars, HD 126587 and CD — 30° 0298, for which 
the uncertainty is as large or even larger than the par¬ 
allax itself. Therefore, the surface gravities of these two 
sta rs were estimated from a c olour-magnitude diagram us¬ 
ing ^hudenBe^^Fa^J^OO^ evolut ionary tracks. As de- 
scribed in ICarcia Perez et Fj ll2005ll . this method was ap¬ 
plied to the entire sample, as a consistency check on the 
log 5 values derived from the parallaxes. A mass of 0.8 Mq 
was assumed for all t he stars and the bo lometric correc¬ 
tions were taken from iBessell et aP lll998ll . 

The metallicity values reported in Table ^ are spec¬ 
troscopic values, and they have been determined from the 
eq uivalent widths of a set o f Fell lines (for more details 
see lCarcia Perez et al ] |2005l) . This same set of iron lines 
was used to constrain the microturbulence, by minimising 
the dependence of the iron abundance on the equivalent 
widths of the lines. 

3.2. Model atmospheres 

In order to model the stellar atmospheres of the pro¬ 
gramme stars, the equation of hydrostatic equilibrium, the 
equation of radiative transfer and the flux-constancy con¬ 
dition, with mixing-length convective flux, were solved si¬ 
multaneously for the set of stellar parameters of each of 
the stars (cf Table [Tt. using an u pdated version of the 
MARCS code llAsohind et al.lll99^ . In order to solve the 
equations (in one dimension and assuming local thermo¬ 
dynamic equilibrium (LTE) conditions), the program re¬ 
quires in addition to the stellar parameters listed above 
also abundances of certain elements such as C, N, O, Ne, 
Na, Mg, etc. These elements need to be taken into ac¬ 
count when computing the number of electrons and opac¬ 
ities in stellar atmosphe res and their relative abu ndances 
were chosen to be solar llOrevesse fc SauvFlll998h . scaled 
to the stellar metallicity; however, the solar oxygen abun¬ 
dance was taken from the ID- LTE analysis of the oxygen 
forbidden line at 630.03 nm ijNisseu et 31.11200^ . and a 
value of 0.4 dex was assumed for the a-enhancement of 
Population II stars. No enhancement was applied to the 
C abundance. 


4. Abundance analysis and uncertainties 

4.1. Lithium abundances 

The Li line at 670.78 nm lies in a spectral region which 
is almost free of other lines, making it possible in all our 
spectra to easily identify a clean wavelength region of ~ 
0.4 nm around the Li line and to unambiguously define 
the stellar continuum. This, together with the fact that 
the blending by a neighbouring Fei line at 670.74 nm is 
insignihcant at our low metallicities, makes the measure¬ 
ment of the equivalent widths very accurate. The widths 


Table 1. Stellar parameters of our programme stars as 
derived bv lGarcia Perez et al.1 ll2005ll 


Star 

Teff 

[K] 

logs 

[cgs] 

[Fe/H] 

[dex] 

^micro 

[km.s“^] 

HD 4306 

4990 

3.04±0.26 

-2.33 

1.5 

HD 26169 

4972 

2.49±0.26 

-2.28 

1.5 

HD 27928 

5044 

2.67±0.43 

-2.14 

1.5 

HD 45282 

5352 

3.15±0.13 

-1.52 

1.1 

HD 108317 

5300 

2.76±0.22 

-2.25 

1.5 

HD 126587 

4712 

1.66±0.50 

-2.87 

1.5 

HD 128279 

5336 

2.95±0.21 

-2.19 

1.5 

HD 200654 

5292 

2.86±0.35 

-2.73 

1.5 

HD 218857 

5015 

2.78±0.36 

-1.79 

1.4 

HD 274939 

5090 

2.79±0.33 

-1.49 

1.3 

BD - 01° 2582 

5072 

2.92±0.40 

-2.11 

1.5 

CD - 24° 1782 

5228 

3.46±0.35 

-2.25 

1.5 

CD - 30° 0298 

5196 

2.93±0.50 

-3.01 

1.5 


were measured by integrating the area under the line us¬ 
ing the IRAF task splat and assuming a Gaussian proHle. 
Direct integration, without the profile assumption, gave 
very similar results. The values are given in Col. 2 of 
Table 121 There is only one star, HD 108317, for which the 
line could not be detected. A value of 0.08 pm was esti¬ 
mated as an upper limit to its equivalent width based on a 
two sigma observational error. The associated abundance 
must therefore be considered an upper limit. 

The lithium line at 670.78 nm is a multiplet but the 
resolution of our observed red spectra is not high enough 
to resolve it completely. The feature was treated as a sin¬ 
gle line when solving the radiative transfer equation. The 
oscillator strength used in our analysis {ioggf = 0.17) 
was the sum of the values for the indivi dual components 
given in Table 3 of ISmith et all lll993ll . ID-LTE abun¬ 
dances were derived from the measured equivalent widths 
using the EQW code from the Uppsala package of stellar 
atmospheres. As a cross-check, we also ran some spectral 
syntheses and found that in general there is a very good 
agreement between the lithium abundances derived by the 
two methods. However, it is well known that the Lii line 
at 670.78 nm forms under non-LTE conditions (NLTE). 
Therefore, NLTE corrections should be applied to the LTE 
abundances. For the stars analysed here, these corrections 
were estimated using the routines and the NLTE results in 
iGsrlsson et af .1 ijl 99411 . The final NLTE Li abundances (i.e. 
LTE Li plus the correction) are listed in Col. 9 of Table|21 
The largest NLTE correction is of the order of 0.20 dex. 
The lithium abundance of HD 108317 and the effective 
temperature of the giant HD 126587 lie outside the range 
of Tes and >l(Li) covered in ICarlsson et alJ l|l994ll . For 
these two stars, NLTE corrections based on their closest 
points in the grid of Carlsson et al. have been adopted. 

The final uncertainty in each lithium abundance was 
calculated by taking into account the uncertainty in the 
equivalent width measurement and in the stellar parame¬ 
ters. According to photon statistics, the associated error 
in equivalent widths for a S/N ~ 500 spectrum is of the 
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Table 2. Lithium abundances of the programme stars. Equivalent widths {Wx) of the Li feature at 670.78 nm are given 
in the second column, followed in the third column by the lithium abundances. Columns 4 to 7 show the sensitivity 
of the lithium abundances to changes in the stellar parameters and in the equivalent width (varied by their typical 
uncertainties). Final errors base d on the square sum o f these values are listed in Col. 8 (a). NLTE Li abundances as 
derived from the multi code o f ICarls son et al ] (l1994^ and the oxygen abundances as derived from the O i lines at 
630.03 nm Icf ICarcia, Perez et a.1 .ll200iT|l are given in the two last columns. 


Star 

IEa 

[pm] 

A(Li) 

LTE 

A±(Li) 

(Ueff) 

AT(Li) 

(logff) 

AT(Li) 

([Fe/H]) 

A±(Li) 

(Wx) 

a 

T(Li) 

NLTE 

±(0) 

LTE 

HD 4306 

1.63 

0.97 

0.11 

0.00 

0.00 

0.02 

0.11 

1.06 

7.14 

HD 26169 

1.44 

0.91 

0.11 

-0.01 

0.00 

0.02 

0.11 

1.08 

7.05 

HD 27928 

1.36 

0.96 

0.10 

-0.02 

-0.01 

0.01 

0.10 

1.12 

7.06 

HD 45282 

1.03 

1.14 

0.10 

0.00 

0.00 

0.02 

0.10 

1.24 

7.72 

HD 108317 

< 0.08 

< 0.00 

0.09 

-0.01 

0.00 

< 0.18 

0.20 

< 0.11 

7.23 

HD 126587 

1.75 

0.76 

0.11 

-0.03 

0.00 

0.01 

0.11 

0.96 

6.32 

HD 128279 

1.06 

1.13 

0.09 

0.00 

0.00 

0.02 

0.09 

1.24 

6.96 

HD 200654 

1.10 

1.09 

0.10 

0.00 

0.01 

0.03 

0.10 

1.20 

< 6.75 

HD 218857 

1.33 

0.90 

0.12 

-0.01 

0.00 

0.02 

0.12 

1.06 

7.33 

HD 274939 

1.44 

1.02 

0.11 

-0.01 

0.00 

0.02 

0.11 

1.17 

7.76 

BD - 01° 2582 

0.63 

0.63 

0.10 

-0.01 

0.00 

0.03 

0.10 

0.77 

7.17 

CD - 24° 1782 

1.12 

1.05 

0.09 

0.00 

0.00 

0.02 

0.09 

1.16 

6.80 

CD - 30° 0298 

1.01 

0.96 

0.09 

0.00 

0.00 

0.02 

0.09 

1.08 

< 6.38 


order of 0.03 — 0.04 pm. Here, in order to take into account 
also the uncertainty due to the continuum placement (al¬ 
though very small) we assumed a total uncertainty of 0.05 
pm in each equivalent width. On the average, this corre¬ 
sponds to abundance errors of the order of 0.02 dex only. 

The sensitivity to the stellar parameters was computed 
by testing the response of the derived lithium abundance 
to changes in TeS (±100 K), in logg (± its la errors which 
are listed in Col. 3 of TableCJ, and in [Fe/H] (± 0.1 dex). 
The influence of these changes on the final Li abundance 
are listed in Cols. 4 to 7 of Table[21 Obviously, the lithium 
abundances are mostly sensitive to the effective tempera¬ 
ture (a 0.1 dex effect). The errors in lithium abundances 
due to errors in other stellar parameters are negligible. 

4.2. Beryllium abundances 

The beryllium abundances are based on synthetic- 
spectrum analysis of the Ben resonant doublet at 313.0-1 
nm. The bluer component (A 313.04 nm) is the stronger of 
the two lines, but it is severely blended (mainly by Cr, V, 
CH and OH lines). The other component of the doublet 
(A313.il nm) is very weak in our programme stars; it lies 
in the wing of a much stronger Ti line but it is otherwise al¬ 
most free of blends. There is no special reason to prefer one 
line over the other. Instead, spectral synthesis of a region 
of 0.3 nm around the two Be lines is mandatory in order 
to properly account for the blends. The ID-LTE synthe¬ 
ses were run with the program bsyn of the Uppsala suite 
of codes. The line list of a tomic lines and CH lines used 
was the same as the one in IPrima s et al.l lll997ll. wherea s 
a new OH line list was used, based on iciilis et alJ ll200lh . 

The elemental abundances were those adopted for the 
model atmospheres, except for those elements with transi- 


Table 3. Our final Be abundances as determined from the 
spectral synthesis fitting of the Be doublet (Col. 2). The 
sensitivity of the result to changes in the stellar param¬ 
eters (corresponding to their uncertainties) are presented 
in Cols. 3 to 5 and their associated a values, calculated as 
in the case of Li, are given in Col. 6. Column 7 lists the 
final NLTE corrected Be abundances. 


Star 

H(Be) 

LTE 

AA(Be)Ayl(Be) AA(Be) 
(Teff) (logs) ([Fe/H]) 

cr 

H(Be) 

NLTE 

HD 4306 

-1.84 

-0.02 

0.17 

0.03 

0.17 

-1.84 

HD 26169 

-1.65 

0.00 

0.15 

0.02 

0.15 

-1.65 

HD 27928 

-1.96 

-0.01 

0.26 

0.02 

0.26 

-2.00 

HD 45282 

-1.15 

0.00 

0.07 

0.03 

0.08 

-1.18 

HD 108317 

< -2.30 

0.03 

0.10 

0.01 

0.10 < -2.26 

HD 126587 

-2.58 

0.02 

0.30 

0.01 

0.30 

-2.37 

HD 128279 

< -2.01 

0.03 

0.11 

0.01 

0.11 < -1.94 

HD 200654 

< -2.06 

0.03 

0.17 

0.01 

0.17< -2.04 

HD 218857 

-1.42 

-0.03 

0.22 

0.03 

0.22 

-1.52 

HD 274939 

-1.20 

-0.03 

0.20 

0.03 

0.20 

-1.28 

BD - 01° 2582 

-1.74 

-0.01 

0.25 

0.03 

0.25 

-1.81 

CD - 24° 1782 

-1.45 

0.00 

0.21 

0.02 

0.21 

-1.54 

CD - 30° 0298 

-2.04 

0.03 

0.27 

0.00 

0.27 

-2.04 


tions in the 0.3 nm wavelength range that was synthesised. 
Once computed, each synthetic spectrum was convolved 
with a profile taking care of macroturbulence, stellar ro¬ 
tation etc., and with a Gaussian with a width correspond¬ 
ing to the resolution of the spectra. The widths of these 
functions were changed together with the stellar abun¬ 
dances until the best fit was achieved. Our final beryllium 
abundances (Col. 2, Table|2Jl correspond to the abundance 
value that best fits both Be ii lines, see Fig. |21 
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Table 4. Logarithmic Li and Be depletion factors for 
the programme stars, based on abundances in metal-poor 
dwarf stars. For Be, factors associated to two different 
Be-vs-0 trends are quoted: tren d in Primas et al . 
in preparation) (Col. 4) an d in B oesgaar d et al.1 
(Col. 5). Predictions based on lPelivannis et al.l l|l99nll ' 
lutionary tracks are given in Col. 3 and Col. 6. 


d. (2005 
n lll999^ 


evo- 


Star 

D(Li) 

obs 

D(Li) 

theor 

D(Be) 

obs 

D(Be) 

Boesgaard 

D(Be) 

theor 

HD 4306 

-1.14 

-1.22 

-1.20 

-0.83 

-0.87 

HD 26169 

-1.12 

-1.22 

-0.91 

-0.51 

-0.87 

HD 27928 

-1.08 

-1.22 

-1.27 

-0.87 

-0.87 

HD 45282 

-0.96 

-1.01 

-1.19 

-1.01 

-0.65 

HD 108317 

< -2.03 

-1.10 

< -1.72 

< -1.38 

-0.75 

HD 126587 

-1.24 

-1.22 

-0.80 

-0.17 

-0.87 

HD 128279 

-0.96 

-1.03 

< -1.09 

< -0.67 

-0.68 

HD 200654 

-1.00 

-1.11 

< -0.95 

< -0.46 

-0.76 

HD 218857 

-1.14 

-1.22 

-1.09 

-0.79 

-0.87 

HD 274939 

-1.03 

-1.22 

-1.34 

-1.17 

-0.87 

BD - 01° 2582 

-1.43 

-1.22 

-1.20 

-0.84 

-0.87 

CD - 24° 1782 

-1.04 

-1.17 

-0.50 

-0.02 

-0.82 

CD - 30° 0298 

-1.12 

-1.19 

-0.54 

0.07 

-0.84 


The formal uncertainty in our Be abundance determi¬ 
nations was derived similarly to what was done in the case 
of lithium. Here we took only into account the sensitivity 
of Be to changes (corresponding to Icr uncertainties) in 
the stellar parameters. We found that the surface gravity 
determination plays a major role in obtaining accurate Be 
abundances, leading to uncertainties H(Be) of about 0.2 
dex. On the other hand, uncertainties associated to effec¬ 
tive temperature and metallicity are insignificant. 

An important source of error in the Be abundance is 
due to the normalisation and continuum location in fit¬ 
ting the synthetic spectra. We estimate that this source 
of error may well contribute errors in A(Be) of about 0.2 
dex. There is also a risk that this error may be systematic, 
e.g. leading to too small Be abundances for the majority 
of the stars. 

The Be abundances given in Table |21 were deter¬ 
mined from the LTE spectral analysis of the Bell UV 
doublet lines at 313.0-1 nm. How ever, according to the 
work of lCfarcfa Lonez et al ] (Il99,^ and Garcfa Perez and 
Kiselman (2005, in preparation), these lines are formed 
in the stellar atmospheres under non LTE conditions. 
Specific NLTE calculations were carried out following the 
procedure described and discussed by Garcfa Perez and 
Kiselman (2005, in preparation), using an adapted version 
of the MULTI code. NLTE Be abundances were determined 
by matching the equivalent width calculated in LTE with 
abundances from the spectrum synthesis with the NLTE 
computed ones. The highest NLTE correction (0.2 dex) 
is found for the giant star, HD 126587. For this star, the 
NLTE Li abundance correction is also among the highest. 


2.5 


2.0 


^ 1-5 
< 

1.0 


0.5 


0.0 

-3.5 -3.0 -2.5 -2.0 -1.5 -1.0 

[Fe/H] 

Fig. 3. Logarithmic lithium abundances in a 12-scale as 
a function of stellar metallicity. Both NLTE (filled circles) 
and LTE (open circles) abundances are shown. Dashed line 
shows a 2.2 constant value as representative of metal-poor 
main-sequence stars in the Spite plateau. 



5. Results and discussion 

Simultaneous knowledge of lithium and beryllium abun¬ 
dances in the same star is a powerful diagnostic of mixing 
activity in stellar interiors. According to standard models, 
halo stars cooler than about 5700 K on the main-sequence 
should show Li depletion (which has already taken place 
during their pre-main-sequence phase). On the contrary, 
metal-poor stars warmer than this temperature all show a 
similar Li abundance, the so-called Sp ite lithium plateau 
I Spite fc SoitelllQ^ iR van et 199911 . According to stel¬ 
lar evolution theory, as stars evolve from dwarfs to giants, 
they get cooler but the surface convection zone deepens 
and reaches layers which were hot enough to “burn” Li. 
The matter poor in Li is mixed with the rest of the con¬ 
vection zone, so that dilution of lithium is expected. 

The surface gravities of our programme stars are typ¬ 
ical of subgiants. Their spectra show weak lithium lines 
which suggest depletion. In one case, HD 108317, the Li 
line was not even detected. Something similar happens to 
the Bell UV doublet lines at 313.0-1 nm in subgiants, i.e. 
they get weaker and more difficult to detect, especially the 
redder line of the doublet which is the cleaner but unfortu¬ 
nately also the weaker of the two. This line (A 313.1 nm) 
was not detected in the observed spectra of HD 108317, 
HD 128279 and HD 200654. When upper limits are given, 
they are constrained by synthetic spectrum fits to both 
lines. 


5.1. Lithium depletion 

As is seen in Figs. 13 and 01 the derived lithium abundances 
of the programme stars, either the LTE (open circles) or 
the NLTE (filled circles) values, are well below the Spite 
plateau value (here taken to be A(Li)NLTE = 2.20, close to 
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Fig. 2. ID-LTE spectral syntheses for a selection of 8 programme stars, in the wavelength range around the Bell UV 
doublet lines at 313.0-1 nm computed for four different Be abundances: the best-fit abundance, abundance changes by 
±0.2 dex around that value, and a very low value of the Be abundance (practically zero). The observed spectra are 
denoted by open circles and the synthetic spectra by solid lines. The dashed line corresponds to the synthetic spectra 
the verv low value of HlBel. 
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Fig. 4. Logarithmic Li depletion factors (top panel) and 
Li abundances (bottom panel) as a function of effective 
temperature. Observed depletion factors and the NLTE 
Li abundances (filled circles) on which these factors are 
based are plotted together with the predicted values (open 
circles). A value of 2.2 was taken as representative of initial 
Li abundance. The triangle denotes the star HD 108317 
with an upper limit on the Li abundance. 
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iBonifacio fc Moliirol l)l997ll l. by about 1 dex. No obvious 
correlation is seen between the lithium abundances and 
metallicity, and the observed difference in 7l(Li) between 
stars at the same metallicity can be fully explained by 
the associated errors, mostly due to uncertainties in the 
effective temperature. 

Assuming that our observed subgiants have left the MS 
with Li abundance values typical of the Spite plateau, we 
can estimate their logarithmic depletion factors using the 
Spite plateau value as the initial value (A(Li)i = 2.2). The 
depletion factors (Col. 2, Table^J based on that MS value 
are compared with dilution predictions (Col. 3, Tabled 
in Fig. 0 as a function of effective temperature, where 
the predi cted values are based on evolutionary tracks pre¬ 
sented in IPeliva.nnis et in ill 99ft for a subgiant star with 
0.775 M 0 and [Fe/H]=-2.3. 

First of all, we note that there is now a trend in the log¬ 
arithmic depletion factor D(Li) (where D{L\) = ^(Li) — 
A(Li)i) with the effective temperature. This is intrinsically 
related to the dre dge-up mechanism and o t her observa¬ 
tional studieSj_e;^JPil 2 £h 2 ^glde_^^ I (|l993|) . Ili bre et all 
l|l999f) and IPe Medeiros et alJ (j2QQC ) have also found a 
similar depletion pattern. If all the stars have left the 
main-sequence with the same Li abundance, then our 
results imply that the cooler subgiants have undergone 
the larger depletions: cooler means more evolved, hence a 
more complete dredge-up mechanism, which is in agree¬ 
ment with standard stellar models (e.g. iDeliva.nnis et a,lj 
Il99f)h . 

Secondly, the observed depletions of our stars seem to 
agree rather well with the predictions, although the ob¬ 
served stars may systematically show somewhat less de¬ 


pletion than the models. We note that BD — 01° 2582 lies 
slightly below the general trend, and that HD 108317, de¬ 
noted by a triangle, seems to be extremely depleted (espe¬ 
cially since its abundance is only an upper limit) in lithium 
compared with the rest of the sample. 

Differences between observed and predicted depletion 
factors seen in Fig.^jare of the order of 0.1 dex, which may 
be insignificant, since the Spite plateau value is known 
with a n accuracy hardly better than 0.1 dex. For in¬ 
stance, 1 I 1999 I 1 suggested a trend of ^(Li) with 

metallicity, and derived a p lateau value of 2.1, whereas 
iMelfedez fc R,amfrej 1 I 2 OO 4 I 1 recently found H(Li)=2.37, 
using a newly derived (high) temperature scale. This is 
just to give an idea of the range of values that one can 
find in the literature, but it should be kept in mind that 
these comparisons are only qualitative, as no attempt has 
been made to put these values on the same scale (i.e. 
taking into account the fact that different authors may 
have used different analytical tools, different model atmo¬ 
spheres, different temperature scales). 

More relevant to our study which deals with sub- 
giant stars is instead the very recent analysis by 
ICharbonnel fc PrimasI ll2005ll . whose sample includes not 
only MS, but also evolved objects. One of their main find¬ 
ings is that the plateau value for their sample of evolved 
stars is, on average, slightly higher (2.235 or 2.259, de¬ 
pending on the lower Tea cut-off defining the plateau re¬ 
gion, 6000 K or 5700 K respectively) than for the dwarfs 
(2.215 or 2.176, for the same Tgff cutoffs, respectively). 
Also, the lithium abundances derived for the evolved stars 
appear to be slightly more scattered than those for the MS 
stars. Although these differences are well within the obser¬ 
vational errors (as quoted by the authors themselves), if 
real they would imply a better agreement between obser¬ 
vations and predictions for our stars (i.e. differences could 
be reduced by up to 0.06 dex, e.g. if we had assumed 
A(Li)i = 2.259 as our initial lithium abundance). 

Clearly, in order to draw firmer conclusions, one should 
in principle compare the observed depletions to more than 
just one set of theoretical models. In practice, this is 
not an easy task, because there have not been simulta- 
neous theoretica l predi ctions for both Li and Be since 
iDelivannis et al.l lll990h . For the purpose of this test, 
Charbonnel {private communication) has kindly provided 
us with three models computed for a similar stellar metal¬ 
licity ([Fe/H] = —2.27) of Deliyannis’ model but for three 
different stellar masses (0.75, 0.80 and 0.85 Mq). As one 
can see from Fig.|S| this new set of models predict different 
depletions for different stellar masses. Also, it is interest¬ 
ing to note that depletion sets in at higher temperatures 
than in the model of Deliyannis: in other words, the same 
amount of depletion is reached at different temperatures in 
the two sets of models. This is not completely surprising, 
as several factors (e.g. initial He content, mixing length, 
opacities, equation of state) may be slightly different in 
these models. 

In order to make as similar compa r ison a s possible to 
the dilution curves of lPelivannis et ahl lll99(ll . we have in- 
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T.„ 

Fig. 5. Observed logarithmic depletion factors (cir¬ 
cles) and predictions based on: evolutionary tracks with 
[Fe/H] = —2.27 and stellar masses of 0.75, 0.80 

and 0.85 Mq (solid lines, from top to bottom) from 
Charbonnel (priv. comm.) and evolutionary track with 
[Fe/H] = —2.30 and stellar mass of 0.775 Mq (dashed 
line) from Deliyannis. For comparison, an interpolated 
model from Charbonnel for a stellar mass of 0.775 Mq 
is also presented (dot-dashed line). 


terpolated Charbonnel’s dilution curves to a stellar mass 
of 0.775 Mq (dot-dashed curve in Fig.^J. Inspection of 
this figure shows that higher discrepancies are found be¬ 
tween observations and predictions when the interpolated 
dilution curve is adopted. It seems that for Charbonnel’s 
dilution curves, masses lower than 0.70 Mq would improve 
the agreement for the majority of our stars. 



[Fe/H] 


-12 



-15 


-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 
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Fig. 6. NLTE and LTE Be abundances as a function 
of metallicity (top panel) and only NLTE values as a 
function of oxygen abundances (bottom panel). The dot¬ 
ted lines correspond to: the Be-vs-Fe (top panel) and 
the Be-vs-0 (bo tt om p anel) correlation for dwarfs from 
iBoesgaard et al.l (Il999ll . The dashed lines correspond 
to the Be-vs-0 correlation for dwarfs from Primas et 
al. (2005, in preparation) (upper line in bottom panel) and 
for our observed subgiants (lower line in bottom panel). 


5.2. Beryllium depletion 


Our beryllium abundances are plotted versus metallicity 
in Fig. El (top panel), with LTE and NLTE values rep¬ 
resented by open and filled circles, respectively. We note 
that, as in the case of lithium, all our stars have beryl¬ 
lium contents much lower than the general Be evolution¬ 
ary trend derived from dwarf star s (shown as a dotted lin e 
of unitary slope in the figure, cf IBoesgaard et ailllOPOll . 
However, since the main process responsible for the for¬ 
mation of beryllium in the Galaxy is cosmic-ray spalla¬ 
tion reactions which involve protons and a particles on 
one side and C and O nucl ei on the other llReeves et alJ 
ll97r)HMenegiiz7,i et al.lll97l|l . it is more instructive to plot 
Be versus oxygen instead of iron. 

The bottom panel in Fig. El shows the correlation be¬ 
tween NLTE beryllium and oxygen abundances. The oxy¬ 
gen abundances w ere derived from the forb idden oxygen 
line at 630.03 nm i Garcfa Perez et 3111200511 . The best fit 
we obtained is a linear fit (dashed line, with slope of 0.8). 
Stars with only upper limits (denoted by triangles) on 
their Be abundances were excluded from the fit, as well as 


CD-24° 1782 (its oxygen abundance comes from OH UV 
lines since its forbidden line was not detected). 

As was also the case for lithium, HD 108317 seems to 
stick out from the rest of the sample, with a beryllium 
content lower than other stars with similar parameters 
and oxygen abundances. The observed pattern could im¬ 
ply that this star has undergone a more severe depletion or 
it may be that its Be abundance before reaching the sub¬ 
giant evolutionary phase was lower than for stars of similar 
oxygen abundances, due to different production scenarios. 
On the contrary, in the case of BD — 01° 2582, the beryl¬ 
lium content is comparable to other stars of the sample 
with similar oxygen abundances, despite its lithium be¬ 
ing lower than our derived trend. This is not remarkable, 
since lithium burns at a lower temperature, i.e. earlier 
than beryllium if the mixing deepens to hotter stellar lay¬ 
ers. 

Furthermore, the stars HD 26169 and HD 27928 have 
very similar O abundances ([O/H] ~ —1-7), but their 
beryllium contents differ by 0.35 dex. We note that this 
difference can be only marginally accounted for by the 
associated observational errors (the final error in our Be 
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Wavelength [nm] 


Fig. 7. The UV Be ii doublet lines in the observed spectra 
of HD 26169 (solid line) and HD 27928 (dot-dashed line). 


abundances is ~ 0.20 dex). As it can be seen from Fig. 0 
where the spectra of the two stars are compared, the dif¬ 
ference in Be abundances may well be real. Having al¬ 
ready measured lithium in both stars, we can now use 
this important information to check if there is any evi¬ 
dence for different mixing history in these two stars. From 
Table[21 we see that HD 26169 and HD 27928 have practi¬ 
cally identical lithium abundances (within 0.05 dex from 
each other) and very similar fundamental stellar param¬ 
eters (cf Table P). One could then expect that they have 
also undergone similar mixing processes: if Be is affected 
by a given transport mechanism, then the Li should have 
changed at least by the same amount. Therefore, the only 
possible explanation, if the Be abundances are proven to 
depart signihcantly in these two stars, is a difference in 
their initial Be abundances: a production more than a des¬ 
truction effect. 

In order to take a look at the global trend of our Be 
abundances, we now compare the Be-vs-0 correlation we 
have derived from our subgiant sample with two linear 
correlations (dotted and upper dashed li ne, bottom panel, 
Fig. 13 derived for dwarfs respectively bv iBoes guard et akl 
(1199911 and more recently also by our group (Primas et 
al. 2005, in preparation). We note that the latter analysis 
finds a non-negligible scatter around the high-metallicity 
end of the Be-vs-0 relation (which overlaps with our more 
metal-rich stars) and that both linear and broken fits are 
under investigation. However, for the purpose of our dis¬ 
cussion, and because the data of Primas et al. are still 
preliminary, we chose to use the simplest of their hts, i.e. 
the linear one. 

Under the assumption that the initial Be abundance of 
all our subgiants was the one found in dwarf stars, we can 
now estimate logarithmic Be depletion factors (I?(Be)), 
based on the NLTE Be abundances reported in Table |21 
It is important to note that because both dwarf correla¬ 
tions reach only [0/H]=—2.0 (whereas our sample reaches 
[0/H]~ —2.5), we decided to extrapolate both hts down 


to our lowest O abundance in order to be able to compute 
the D(Be) factors for the entire subgiant sample. The fac¬ 
tors thus derived suggest that our stars have depleted Be 
by approximately 1 dex, as they have evolved off of the 
main sequence. 

Similarly to Li, we compare our estimates (Col. 4, 
Table EJ with predictions (Col. 6, Tabled based on de¬ 
pletion calculations kindly provided by C. Deliyannis (pri¬ 
vate communication). As it is shown in the upper panel 
of Fig. H the standard models predict a Be dilution in 
subgiants which is expected to depend on effective tem¬ 
perature. However, our depletion factors do not show any 
obvious sign of such a dependence, although the observa¬ 
tional or cosmic scatter may be too large to disclose such 
a tendency. 

Because of the different slopes in the dwarf correla¬ 
tions we have used as a comparison, the emerging trends 
are quite different. Be seems to be systematically more 
depleted in our subgiant sample (on average by 0.2 dex) 
than expected from the depletion calculations (compare 
filled with empty circles) if the Primas et al. corr elation is 
used. Instead, when the iBoesgaard et al ] (Il999h relation 
is considered, the difference between observed and pre¬ 
dicted depletion values is reduced (Cols. 5 and 6, Tabled, 
and for some objects we even observe a reversed behavior, 
i.e. the models seem to overproduce the amount of deple¬ 
tion that Be has already undergone. Note, however, that 
the relation of iBoesgaard et al~l l|l999l) is based on oxygen 
abundances derived from OH UV lines instead of the for¬ 
bidden oxygen line at 630.03 nm. For these reasons and 
because we and Primas et al. have used similar methods 
to derive abundances and stellar parameters, we believe 
that it is more appropriate to carry out the comparison 
between MS and subgiant stars in terms of the Primas et 
al. dwarf relation. 

Furthermore, the derived depletion factors do not fol¬ 
low a clean trend but they are quite scattered, especially 
in a range of 150 K around 5050 K, where the spread in the 
Be depletion factors is considerable (with a maximum dif¬ 
ference of the order of 0.4 dex). Although this seems to be 
at odds with the constant depletion value expected from 
dilution at these temperatures, the observational errors as¬ 
sociated to our Be abundances are quite large. Therefore, 
the observed scatter cannot be taken as a clear indication 
of physical differences between the stars. 

Possibly, the systematic discrepancies between ob¬ 
served and predicted depletion factors for Be may be the 
product of our assumption that initial Be abundances are 
those found in dwarf stars. We explore this possibility by 
comparing estimates with predictions in a I?(Be)-vs-[0/H] 
plot (bottom panel in Fig. 0. We do not see, however, sig- 
nihcant differences in depletion discrepancies for different 
[O/H]. This suggests that if the assumed values for the ini¬ 
tial Be abundance were wrong, the relative errors would 
be independent of [0/H]. 

According to the tracks of Charbonnel and as it was 
the case for Li, Be depletion depends on the stellar mass. It 
is also the case here that the Charbonnel’s dilution curve 
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Fig. 8. Logarithmic Be depletion factors (top panels) and 
Be abundances (bottom panels) as a function of effective 
temperature (left panels) and oxygen abundance (right 
panels). Our estimates (filled circles) are compared with 
the predictions (open circles). Stars with upper limits on 
Be abundance are denoted by triangles. 


Fig. 9. The differences between logarithmic depletion 
factors for Li and Be as a function of temperature (top 
panel) and oxygen abundance (bottom panel). Estimates 
(filled circles! are com pared with predictions based on 
IPelivannis et al.l lIlflQflll (open circles). 


interpolated for a stellar mass of 0.775 Mq produces more 
depletion than the dilution curve of Deliyannis (see right 
panel in Fig|3. This clearly helps to reduce the differ¬ 
ences between observations and predictions, although this 
is not enough to eliminate them completely. For most stars 
in our sample and unlike the case of Li, it would be nec¬ 
essary to resort to dilution curves corresponding to stellar 
masses higher than 0.775 Mq, in order to reconcile our re¬ 
sults with the predictions. Furthermore, we note that dif¬ 
ferences between observations and predictions are much 
higher than the differences between the dilution curves 
shown here. 

5.3. Comparison between Li and Be depletions 

Lithium is the most fragile of the three light elements. If 
the Be content of a star is depleted, it is expected to have 
already depleted all, or a great fraction of its lithium. 

According to standard model predictions (open circles, 
Fig. El), Li depletion should be around 0.3-0.4 dex higher 
than Be depletion. Despite this is a feature common to 
both sets of dilution curves we have used in our compar¬ 
isons, most of our observed subgiants (more than 50%) are 
characterised by much larger Be depletions (compared to 
their amount of Li depletion). The stars with a detected 


Li line but with only upper limits for the Be abundance 
(HD 128279 and HD 200654) are of special interest in this 
respect. 

As already discussed in Sect. 5.1, the agreement be¬ 
tween observed and predicted D(Li) values is at the level 
of 0.1-0.15 dex, and it holds for both sets of dilution 
curves (i.e. Deliyannis and Charbonnel) in the case of the 
cooler objects of our sample. Thus, it must be beryllium 
that sticks out: models predict a much smaller amount 
of depletion, between 0.25 and 0.65 dex, depending on 
the star under consideration and also on the set of mod¬ 
els (for Be, Charbonnel’s dilution curves predict a larger 
amount of depletion for the warmer stars of our sample). 
Furthermore, one feature that could be worth further in¬ 
vestigating is the metallicity dependence of the theoretical 
dilution curves: in fact, we noted that the best agreement 
between observed and predicted depletion factors (when 
both Li and Be are taken into account simultaneously) is 
obtained for the two most metal-poor stars of our sam¬ 
ple (CD — 24° 1782 and CD — 30° 0298), while the worse 
agreement is obtained for the most metal-rich stars of our 
sample (HD 45282 and HD 274939). 

Of course, another possibility is that our Be abun¬ 
dances are systematically too low by about 0.2 dex as 
a consequence of a too low spectrum normalisation (con¬ 
tinuum location). This may be directly seen in Fig. El 
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where we would not need to change the continuum level 
very much in order to get a good fit of the observed Be ii 
lines with the lower solid lines. With such Be abundance 
uncertainty in mind, we should refrain from any defini¬ 
tive statement on departures from the standard depletion 
calculations, as well as on the possibility that the Be de¬ 
pletion could be larger than that of Li. However, we can 
make some qualitative speculations. 

If our abundances are taken at face value (i.e. without 
any systematic errors in the Be abundances), we could 
still reconcile the Li and Be depletion data by advocating 
a higher Li depletion, e.g. assuming a higher initial lithium 
of ^(Li) = 2.6. As a matter of fact, the primordial lithium 
abundance infe rred by the baryonic density obtained from 
WMAP data dSpergel et al.l Efio^ in combin ation with 
stand ard Big Bang nucleo-synthesis (SBBN, ICoc et aP 
I2nn4l) is ^(Li) = 2.6. However, since much lower values 
are consistently observed in metal-poor dwarf stars (and 
most of them cluster around H.(Li) = 2.2), such a scenario 
is unlikely. Mixing models like those includi ng rotation- 
ally induced mixing llPinsonneault et alJ2002^ indeed sug¬ 
gest some degree of depletion in these stars with respect 
to the primordial value, the magnitude of which strongly 
depends on stellar properties such as rotational veloci¬ 
ties, masses, temperatures, metallicities etc. However, the 
absence of any significant scatter in the Li abundances 
measu red in warm halo stars does not favour suc h a sce¬ 
nario ( Bonifacio fc Mola^ll997HR,van et al.lll999fl . and it 
is not obvious whether mixing processes are able to pro¬ 
duce such high Li and Be depletions. We note that theore¬ 
tical computations including rotation mixing and gravity 
wave s seem to be able to deplete Li in a very homogeneous 
way l|Talon fc Charbonnelll2004 . as required by what is 
observed in warm halo dwarf stars. Unfortunately, Be de¬ 
pletions have not yet been derived and inspected. The 
luminosities of our metal-poor subgiants suggest that sev¬ 
eral of them are starting to ascend the red giant branch 
so they may already have suffered a first dredge-up. If 
this process is not well described by standard theory, the 
predictions of both Li and Be may be underestimated. 


6. Conclusions 

As stars evolve off the main sequence, their surface convec¬ 
tion zone starts to move inwards reaching first Li depleted 
and later Be depleted matter. The mixing between these 
outer and inner layers reduces the surface abundances of 
these light elements. Subgiants are expected to show signs 
of such depletion in their Li and Be abundances as com¬ 
pared with MS values. Our derived Li abundances lie sig¬ 
nificantly below the Spite lithium plateau. The depletion 
is approximately one order of magnitude. A trend with ef¬ 
fective temperature may be traced: the warmer stars show 
the higher Li abundances. In general, these abundances 
agree well wit h the standard p r edicti ons from evolution¬ 
ary models of iDelivaunis et ^ ()l 99flll with the provision 
that the observations may show somewhat less Li deple¬ 
tion than expected. These differences can be reduced to 


zero if a value for the Spite plateau 0.1 dex higher were 
assumed instead, which is quite possible given the uncer¬ 
tainty of its value. 

As far as Be is concerned, observations of subgiants 
show that abundances of Be and O are correlated, as in 
the case of dwarf stars. Be depletion estimates depend on 
the assumed Be-vs-0 trend of dwarfs, e.g. in our case, 
that of Boesgaard et al. (1999) and of Primas et al. (2005, 
in preparation). However, since similar stellar parameters 
and analytical methods were used here and by Primas et 
ah, we have considered the latter source more appropri¬ 
ate for estimating our Be depletion factors. What emerges 
is that observed and predicted depletions in our subgiant 
sample are comparable but in some conflict; the predicted 
Be depletion is typically 0.3 - 0.4 dex smaller than what we 
observe. We note that a higher initial amount of Li would 
help reconciling our derived D{L\) and D(Be) depletion 
factors, so that D(Li) is larger than D(Be) for all our stars. 
This is challenged by the robustness of the Spite’s plateau 
value of A(Li) ~ 2.2. A wrong temperature scale could 
affect this value, but not enough to reconcile our II(Li) 
and II(Be) values. This would require ^(Li) ~ 2.6 which 
is an unlikely value for the initial Li content of dwarf stars 
but surprisingly matches the WMAP-I-SBBN predictions 
for the primordial Li content. An alternative would be de¬ 
pletion mechanisms acting more efficiently for Be than for 
Li, mechanisms that are not easy to envisage. However, 
we note that 0.3 dex could also be a (conservative) esti¬ 
mate of the uncertainty in our Be abundances due to the 
continuum placement. Hence, should we have systemati¬ 
cally performed our spectrum normalisation too low, then 
a better agreement between standard predictions and ob¬ 
servations is obviously obtained. 

Similar conclusions can be derived when Charbonnel’s 
dilution curves are used: namely relatively good agree¬ 
ment between observations and predictions of Li depletion, 
but somewhat worse in the case of beryllium. According 
to these dilution curves, our results on Li suggest that 
the observed subgiants should have masses smaller than 
0.75 Mq, whereas our observed Be depletion factors sug¬ 
gest that a better agreement could be achieved for stellar 
masses larger than 0.85 Mq. 

Further disentangling between light element observed 
abundances and theoretical predictions require the avail¬ 
ability of a grid of dilution curves, covering the whole 
range of metallicities of our sample and different masses 
(since there seem to be a dependence on the stellar mass), 
as well as higher quality data, in order to reduce the error- 
bars associated to individual abundance measurements. 
Should our Be abundances and their depletion factors be 
confirmed, then this will have important implications on 
the modelling of the interiors of these stars. 
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